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We investigate the enhancement of the Kondo effect in quantum dots with an even number of 
electrons, using a scaling method and a mean field theory. We evaluate the Kondo temperature 
Tk as a function of the energy difference between spin-singlet and triplet states in the dot, A, and 
the Zeeman splitting, Ez- If the Zeeman splitting is small, Ez <C Tk, the competition between the 
singlet and triplet states enhances the Kondo effect. Tk reaches its maximum around A = and 
decreases with A obeying a power law. If the Zeeman splitting is strong, Ez 2> Tk, the Kondo 
effect originates from the degeneracy between the singlet state and one of the components of the 
triplet state at — A ~ Ez- We show that Tk exhibits another power- law dependence on Ez- The 
tIh . mean field theory provides a unified picture to illustrate the crossover between these regimes. The 

enhancement of the Kondo effect can be understood in terms of the overlap between the Kondo 
' resonant states created around the Fermi level. These resonant states provide the unitary limit of 

the conductance G ~ 2e 2 jh. 
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The Kondo effect observed in semiconductor quantum dots has attracted a lot of interest IMi In a quantum dot, 
the number of electrons N is fixed by the Coulomb blockade to integer values and can be tuned by the gate voltage. 
Usually the discrete spin-degenerate levels in the quantum dot are consecutively occupied, and the total spin is zero 
J> ■ or 1/2 for an even and odd number of electrons, respectively. The Kondo effect takes place only in the latter case. 
The spin 1/2 in the dot is coupled to the Fermi sea in external leads through tunnel barriers, which results in the 
formation of the Kondo resonant state at the Fermi leveLETQ The conductance-through the dot is enhanced to a value 
of the order of e 2 /h at low temperatures of T <C Tk (Kondo temperature) .lTEj This is called unitary limit. When N 
, is even, there is no localized spin and thus the Kondo effect is not relevant. 
__, ' Recently Sasaki et al. has found a large Kondo effect in so-called "vertical" quantum dots with an even N$3 
The spacing of discrete levels in such dots is comparable with the strength of electron-electron Coulomb interaction. 
■ Hence the electronic states deviate from the simple picture mentioned above. lj£3 If two electrons are put into nearly 

S degenerate levels, the exchange interaction favors a spin triplet (Hund's rule).t3 This state is changed to a spin singlet 
by applying a magnetic field which increases the level spacing. Hence the energy difference between the singlet and 
triplet states, A, can be controlled experimentally by the magnetic field. The Kondo effect is significantly enhanced 
around the degeneracy point between the triplet and singlet states, A = 0. Tuning of the energy difference between 
the spin states is hardly possible in traditional Kondo systems of dilute magnetic impurities in metal and thus this 
£j , situation is quite unique to the quantum dot systems. ._. _ . 

The Kondo effect in multilevel quantum dots has been investigated theoretically by several g:pups.EZn 2 i3 They have 
shown that the contribution from multilevels enhances the Kondo effect. In our previous paper J 2 -!! we have considered 
the experimental situation by Sasaki et al. in which the spin-singlet and triplet states are almost [degenerate. We have 
calculated the Kondo temperature Tk as a function of A, using the poor man's scaling method.E3~o We have shown 
that Tk(A) is maximal around A = and decreases with increasing A obeying a power law, Tk(A) oc 1/A 7 . The 
exponent 7 is not universal but depends on a ratio of the initial coupling constants. Our results indicate that the 
Kondo effect is enhanced by the competition between singlet and triplet states, in agreement with the experimental 
findings S3 

We have disregarded the Zeeman splitting of the spin-triplet state, —EzM (M = 0,r±l is z-component of the total 
spin S = 1), since this is a small energy scale in the experimental situation, Ez <C TkJii 
et al. have studied another situation where the Zeeman effect is relevant, Ez 3> Tk- 



J^ustilnik et al. and Giuliano 
Usually the Zeeman effect 



lifts off the degeneracy of the spin states and, as a result, breaks the Kondo effect. They have found that the Kondo 
effect can arise from extra degeneracy between one of the components of the spin-triplet state, \SM) = |11), and a 
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singlet state, |00), if the value of A is tuned to fulfill that Ez = —A. Their mechanism might explain some other 
experimental results of the Kondo effect in quantum dots under high magnetic fields D'EH 

The purpose of the present paper is to construct a general theory for the enhancement of the Kondo effect in 
quantum dots with an even number of electrons, for various values of A and Ez- We adopt the poor man's scaling 
method along with the mean field theory. It is well known that the characteristic energy scale of the Kondo physics, 
the Kondo temperature Tk, is determined by all the energies from Tk up to the upper pjtoff.Ljtl By the scaling method, 
we can evaluate Tk (its exponential part at least) by taking all the energies properly.E5l2il When Ez is negligible, the 
energies from A to the upper cutoff would feel fourfold degeneracy of the dot states, |1-S*0 [M = 0, ±1) and |00), 
which enhances the Kondo temperature. With increasing A, Tk decreases by a power la^pp We extend our previous 
calculations to the case of Ez = -A > Tk which has been discussed by Pustilnik et alEB and Giuliano el a/.t3 We 
take into account the energies not only from Tk to Ez, where only two degenerate states |11) and |00) are relevant, 
but also from Ez to the upper cutoff, where the dot states seem fourfold degenerate. The latter energy region has 
been neglected in Refs. ^5] and |2^. In consequence we find a power law dependence oi-Tk on Ez again. 

The mean field theaty of the Kondo effect was pioneered by Yoshimori and SakuraicSl and is commonly used for the 
Kondo lattice model.E3 It is useful to capture main qualitative features of the Kondo effect; renormalizability at the 
scale of Tk, resonances at the Fermi level, and resonant transmission. The simplicity and universality of the mean field 
theory have driven us to apply it to the problem in question. Generally the Kondo effect gives rise to a many-body 
ground state which consists of the dot states \SM) = fg M \0) and the conduction electrons IIc^lO). The total spin of 
this ground state is less than the original spin S localized in the dot. The binding energy is of the order of the Kondo 
temperature Tk- We take into account the spin jCouplings between the dot states and conduction electrons, (fg M Ck a ), 
by the mean field, neglecting their fluctuations .E2I These spin couplings give rise to resonant states around the Fermi 
level n with the width of the order of Tk- The conduction electrons can be transported through the resonant levels, 
which yields the unitary limit of the conductance G ~ 2e 2 /h. For our study, the mean field calculations have the 
following advantages, (i) The enhancement of Tk by the competition between the singlet and triplet states can be 
directly understood in terms of the overlap between their Kondo resonant states, (ii) The power law dependence of Tk 
on A or Ez is obtained, which is in accordance with the calculated results by the scaling method, (iii) The mean field 
calculations are applicable to the intermediate regions where two of Tk, A, and Ez, are of the same order. The poor 
man's scaling method hardly gives any results in these regions. Hence we can examine the whole parameter region of 
A and EsJay the mean field theory. The disadvantage of the mean field calculations is that they only give qualitative 
answers.E2l Hence the mean field theory and scaling method are complementary to each other for understanding the 
Kondo effect. on 

We shall discuss the relation of our approach to the renormalization group analysis of the multilevel Kondo cffcctoo 
Our model effectively reduces to the one with two channels in the leads and spin-triplet (and singlet) state in the dot. 
The ground state of this model would be believed to be a spin singlet, which corresponds to the full screening of the 
dot spin. The poor man's scaling approach and our mean field theory, however, show a tendency to the formation 
of the underscreened Kondo ground state with spin 1/2. We should mention that the exact ground state can not be 
determined within the limits of the applicability of these-approaches. Pustilnik and Glazman have recently proposed 
a different model for the "triplet-singlet Kondo effect. "Eil In our notations, they set C\ = \/2, C'2 = in Eq. (1 
the singlet state. Their model can be directly mapped onto a special case of the two-impurity Kondo model, 
which the ground state is a spin singlet. We are concerned about the case of C\ « C2, and we find that the difference 
between C\ and C2 reduces as a result of the renormalization.El This suggests that the case considered in Ref. ||l] is 
by no means a generic one. 

This paper is organized as follows. Our model is presented in the next section. In section III, we rederive Xk(A) 
when_the Zeeman splitting is irrelevant, using the poor man's scaling method, in a simpler form than our previous 
work.B Then we extend our calculations to the case of Ez = — A 3> Tk. The section IV is devoted to the mean field 
theory for the Kondo effect in quantum dots. First we explain this theory for the usual Kondo effect in a quantum 
dot with S — 1/2. Then we apply the mean field scheme to our model with an even number of electrons in the dot. 
The conclusions and discussion are given in the last section. 



II. MODEL 



We are interested in the competition between the spin-singlet and triplet states in a quantum dot. To model the 
situation, it is sufficient to consider two extra electrons in a quantum dot at the background of a singlet state of all 
other N — 2 electrons, which we will regard as the vacuum |0). These two extra electrons occupy two levels of different 
orbital symmetry.E3 The energies of the levels are E\ and £2- Possible two-electron states are (i) the spin-triplet state, 
(ii) the spin-singlet state of the same orbital symmetry as the triplet state, 1/V2(c4|<4 1 — ^ii^2t)I^)> ana - ( m ) ^ w0 
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singlet states of different orbital symmetry, d^dj, |0), d^dj ,|0). Among the singlet states, we only consider a state 
of the lowest energy, which belongs to the group (iii). Thus we restrict our attention to four states, \SM): 

iii) = d t 1T 4 T | ) (i) 

|10) = -^(dl T dJ ; + dl ; dt T )|0) (2) 
|l-l) = d t u d|jO) (3) 
|00) = -^(dd^d^ - Cad+^IO), (4) 

where dt creates an electron with spin a in level i. The coefficients in the singlet state, C\, Ci (|Ci| 2 + | C 2 1 2 = 2), 
are determined by the electron-electron interaction and one-electron level spacing 6 = £2 — £i< We set C\ = C2 = 1. 
This is the case for 5 = O.tfl Although C\ ^ C2 in general, the scaling, analysis shows that the Kondo temperature is 
the same as that in the case of C\ = C% = 1, apart from a prefactorO The energies of the triplet state are given by 

E s =i,m = E s =i - E Z M (5) 

and the energy of the singlet state is denoted by E 00 . We define A by 

A = E 00 -E S=1 . (6) 

The energy diagram for the spin states is indicated in Fig. 1(a). 

The dot is connected to two external leads L, R with free electrons being described by 

-n loads — ^ Z—/ k c a,kcr t - a ,ka^ 
a—L : R kai 

where c£'l ( c a\a) ^ s ^ ne creation (annihilation) operator of an electron in lead a with momentum k, spin a, and 
orbital symmetry i (= 1, 2). The density of states v in the leads remains constant in the energy band of [— D, D}. The 
tunneling between the dot and the leads is written as 

H *= E E(^.4 4 iU-+H.c)- 

a—L.R kai 

We assume that the orbital symmetry is conserved in the tunneling processes!!] To avoid the complication due to the 
fact that there are two leads a = L, R, we perform a unitary transformation for electron modes in the leads along the 

lines of Ref. |; c£> = [V£^ M + V^c^/Vi, cf = {-V R>i cf M + VlJ^/K, with V t = y/\V L ,i\* + \V Rii \*. The 
— H) 

modes & ka are not coupled to the quantum dot and shall be disregarded hereafter. Then i?i oa ds and Ht are rewritten 



as 

.(0„Wt„(0 

r 

kai 



^lea ds = E4 i) 4f4^ (7) 

^T = E^(42 td -+ H - c -)- (8) 

kai 

We assume that the state of the dot with N electrons is stable, so that addition/extraction energies, E ± = 
E(N ± 1) — E(N) =F n where /1 is the Fermi energy in the leads, are positive. We are interested in the case where 
E ± ^> |A|, 5 and also exceed the level broadening T l = irvV? and temperature T (Coulomb blockade region). In 
this case we can integrate out the states with one or three extra electrons. This is equivalent to Schrieffer- Wolff 
transformation which is used to obtain the conventional Kondo model JD'B We obtain the following effective low-energy 
Hamiltonian 

H e S — -ffloads + -^dot + + iJ 5 - 1 ^ -f H' eS . (9) 

The Hamiltonian of the dot -f/dot reads 

Hdot = E EsAifsAffsM, (10) 
S,M 
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using pseudo-fermion operators fl M (/sm) which create (annihilate) the state \SM). The condition of 



SM 



(11) 



SM 



should be fulfilled. The third term if 5-1 represents the spin flip processes among three components of the spin-triplet 
state. This resembles a conventional Kondo Hamiltonian for S = 1 in terms of the spin operator S 



H s=i - V V rd^ r «t r W , o J0t r W , 9 f„(0U<) _ r Mt (*K 



fcfc' i=l,2 



£ £ J« [^(AVxo + AVi-Oc^/cg + V2(fl fn + / 1 t _ 1 /io)c 



(0t„(0 

fe'T c fcl 



fcfc' i=l,2 



+(/i t i/n-/ l t -i/i-i)(4v T 1 



c fcT c fc'i c fci./ 



(12) 



The exchange coupling jW is accompanied by the scattering of conduction electrons of channel i. The fourth term 
jjS=i^o m describes the conversion between the spin-triplet and singlet states accompanied by the interchannel 
scattering of conduction electrons 



H S^o = J- j[V2(flj 00 fUi-i)c™c§ + 72(4/11 - fljoo)^, 

fcfc' 

-Uio/oo + Joo/ioJ(,c fc , T c feT - c fc , x c fci ) + (1 

The coupling constants are given by 



)t„(2) 



T/2 

~2£ c ' 

^ 2 



J = 



2E r 



(13) 

(14) 
(15) 



where 1/E C = 1/E + + 1/E~ . Note that J 2 = J^'J&K The last term iJ^ ff represents the scattering processes of 
conduction electrons without any change of the dot state and is not relevant for the current discussion. The spin-flip 
processes included in our model are shown in Fig. 1(b). 



III. SCALING CALCULATIONS 

In this section we calculate the Kondo temperature Tk using the poor man's scaling technique .E1~§1 By changing 
the energy scale (bandwidth of the conduction electrons) from D to D — \dD\, we obtain the scaling equations using 
the second-order perturbation calculations with respect to the exchange couplings, J' 1 ', </' 2 ', and J. With decreasing 
D, the exchange couplings are renormalized. The Kondo temperature is determined as the energy scale at which the 
exchange couplings become so large that the perturbation breaks down. 



A. In the absence of Zeeman effect 



When the Zeeman splitting is. small and irrelevant, E% -C Tk, we obtain a closed form of the scaling equations for 
jW, j( 2 ), and J in two limits.Eil (i) When the energy scale D is much larger than the energy difference |A|, H^ ot can 
be safely disregarded in H c g. The scaling equations can be written as 

J (1) J \ „ f J {1] J 



= - 2 M ~r r(2) ■ (16) 



dlnD V J J {2) )~ \ J J { 

(ii) For fl« A, the ground state of the dot is a spin triplet and the singlet state can be disregarded. Then and 
j( 2 ) evolve independently 
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d\nD 



(17) 



whereas J does not change. 

In the case of |A| <§; Tk, the scaling equations (|l6|) remain valid till the scaling ends. The matrix in Eq. Jl6| ) has 
eigenvalues of 



J± = (J (1) + J (2) )/2± y / (j( 1 ) - J( 2 )) 2 /4 + J 2 

= jw + jw 0. 



(18) 



The larger one, J+, diverges upon decreasing the bandwidth _D and determines Tk 

T K (0) = A,exp[-l/2z/J + ] 

= A,exp[-l/2*/(J (1) + J (2) )]. 

Here -Do is the initial bandwidth, which is given by VE+E^E 

When A > Dq, the scaling equations ([17]) work in the whole scaling region. This yields 

T K (oo) = T> cxp[-l/2iyJ (1) ] 

for > j( 2 \ This is the Kondo temperature for spin-triplet localized spins.El 

In the intermediate region of Tk(0) <C A -C Dq, the exchange couplings develop by Eq. (|T(| ) for D 3> A. Around 
T> = A, J saturates while J*- 1 - 1 and J*- 2 - 1 continue to grow with decreasing D, following Eq. (|l7|) for D -C A. We match 



(19) 



(20) 



the solutions of these scaling equations at D ~ A and obtain a power law of Tk(A) 

T K (A)=TK(0)-(TK(0)/A) tan29 : 

with 

tan (9 



(21) 



J/[ V /(J( 1 ) - J( 2 )) 2 /4 + J 2 + (J (1) - J (2) )/2] 



(22) 



for jW > J^ 2 ). Here (cosfl, sin#) T is the eigenfunction of the matrix in Eq. (|16[) corresponding to J+. 6> ~ for 
j(i) > j(2) anc j = f or j(i) _ j(2)_ In general, < 6 < tt/4 and thus < tap 6> < 1. 

Finally, for A < 0, all the coupling constants saturate and no Kondo effect is expected, provided |A| ^> Tk(0). 
Thus Tk quickly decreases to zero at A ~ — Tk(0). The Kondo temperature as a function of A is schematically shown 
in Fig. 2(a). 



B. Case of E z = -A 

When Ez = — A-qfeiie energies of states |00) and |11) are degenerate. Then the Kondo effect is expected even 
when |A| ^> Tk(0)£3 : e3 In this subsection we evaluate Tk in this special case of Ez = —A by the poor man's scaling 
method. 

(i) For the energy scale of D 3> |A| = Ez, -ffdot can be disregarded in i? c ff- The exchange couplings, JW, j( 2 \ 
and J, evolve following Eq. (|l6|). (ii) In another limit of D <C |A| = Ez, only the states 1 00) and |11) are relevant. In 

H c fi, 

rr|00),|ll> _ V- [1 , ,f , \/-(0t„(0 JO t JO \ , 1 I f t f , Af uJOtJO JOtJ0\' 

-"eff - 2^ [2 * i wii/11 - Joo/ooj( c fc , T c feT - c fe , x c fci j + - Jl )'{} n Ji 1 + JooJao){c k ^c k ^ -c k>i c H ) 



kk' i=l,2 



E ^[/iVooc^Mf + /oVn^i^i? + (1 - 2)1 . (23) 



Js^ = Jc^ = J^- 1 initially. The scaling procedure yields 



r(0 
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^ J s (l) = -Ah>J 2 



iS^-^+^v, (24) 

and J^' fdo not change. These scaling equations are nearly equivalent to those of the anisotropic Kondo model with 
S = l/2,Has pointed out in Refs. ||,[|| 

When |A| = Ez > D Q , the scaling equations ( |24| ) remain valid in the whole scaling region. This yields the Kondo 
temperature 

T K (oo) = D o exp[-A(0)/2z/(jW + j( 2 ))] (25) 

with 



Jll*(£x) (0<e<n/8) 
\ | tan" 1 A (tt/8 < 9 < vr/4), 



A{6) = { * m \?=X) ^<^«/°; (26) 



where A = yj cos 40|. A(0) decreases monotonically with increasing 9. A(9) — ► oo as — > 0. A(7r/8) = 2 and 
,4(71-/4) = tt/2. When jW + J< 2 ) is fixed, T K (oo) is the largest at JW = j( 2 ) (9 = 7r/4) and becomes smaller with 
decreasing j( 2 )/jW(= tan 2 9). 

In the intermediate region, Tk(0) <C I A| = Ez -C A), we match the solutions of Eqs. ( |l6|) and ( ^) at £> ~ |A|. We 
obtain a power law 

T K (A) = T K (0) • (TkW/IAI)^- 1 . (27) 
Figure 2(b) shows the behaviors of 7k(A) in the case of Ez = —A. 



IV. MEAN FIELD CALCULATIONS 

A. Kondo resonance for spin 5=1/2 

To illustrate the mean field theory for the Kondo effect in quantum dots, we begin with the usual case of S = 1/2. 
We assume that one level (Eg) in a quantum dot is occupied by an electron with spin either up or down (a =T,J.)- 
The effective low-energy Hamiltonian is 

h = E £ * c L c ^ +Y, E °tif° + J EE flf^L^, (28) 

ka o kk' <y,a' 

with the constraint of 

/f/T + f\fi = 1- (29) 



For electrons in leads L, R, we have performecLa unitary transformation of c k(J = (V^c i fc(T + V^c R ^ ka )/ \J\Vl\ 2 + \Vr\ 2 
where V a is the tunneling coupling to lead aB The last term in Eq. ( pj| ) represents the exchange coupling between 
S = 1/2 in the dot and conduction electrons (Appendix A). 
In the mean field theory, we introduce the order parameter 

< s ) = 7fE«A tc fcT> + (/N» (30) 

k 

to describe the spin couplings between the dot states and conduction electrons. The mean field Hamiltonian reads 

h mf = E e * c Lc feCT +j2 E °f f J° - E(^ J ( s > c L/- + h.c.) + 2j|(s)i 2 +x(y i flu - 1 ) . (31) 



The constraint, Eq. Q29|), is taken into account by the last term with a Lagrange multiplier A. By minimizing the 
expectation value of Hmf, (S) is determined self-consistently (see Appendix A). 



G 



In the absence of the Zeeman effect, = E± = Eq. The mean field Hamiltonian, Hmf, represents a resonant 
tunneling through an "energy level," E = E + A, with "tunneling coupling," V = —y/2J(S). V provides a finite 
width of the resonance, Ao = m/\V | 2 , with v being the density of states in the leads. The constraint, Eq. (p9|), requires 
that the states for the pseudo-fermions are half-filled, that is, E = fi. Hence the Kondo resonant state appears just 
at the Fermi level fi, as indicated in the inset (A) in Fig. 3(a). The self-consistent calculations give us the resonant 
width 



A = ttv V2J(E) = D Q exp[-l/2vJ}. 



(32) 



This is identical to the Kondo temperature 7k. 

In the presence of the Zeeman splitting, E^ = Eq — Ez and E± 
spin-up and down electrons, Ey± = Ey^ + A. The constraint, Eq. (| 
The resonant width A is determined as 



Eq + Ez ■ Hence the resonant level is split for 
f), yields Eq + A = \i (see inset (B) in Fig. 3(a)). 



A 



(33) 



where Ao is given by Eq. (|32|). The Kondo temperature is evaluated by this width, 7k (-Ez) = A. 2k decreases with 
increasing Ez and disappears at Ez = 7k (0), as shown in Fig. 3(a). 

The conductance G through the dot is expressed, using r Q = iri/\V a \ 2 , as 



G 



2e 2 4r L r fl 
h (r L + r R y 



\Tk(Q)) 



(34) 



This is the conductance in the unitary limit for Ez = 0. Figure 3(b) presents the Ez dependence of the conductance. 
With increasing Ez, the splitting between the resonant levels for spin up and down becomes larger. In consequence 
the amplitude of the Kondo resonance decreases at /i, which reduces the conductance. 



B. Kondo resonance in the present model 



Now we apply the mean field theory to our model which has the spin-triplet and singlet states in a quantum dot. 
The spin states of the coupling to a conduction electron are (S = 1) <g> (S = 1/2) = (S = 3/2) ® (S — 1/2) for the 
former, and (S = 0) ® (S = 1/2) = (S = 1/2) for the latter (Appendix B). To represent the competition between the 
triplet and single states, therefore, the order parameter should be a spinor of S = 1/2. It is (£) where 



2 = £ 



cos ip 



,S ip (V^cff + /iVil) + Sin P/OO 



'k.l + flO C k1 



I ^3 - sin pfl 



,(2) 



(2) 



(35) 



for > j( 2 \ A mode of the largest coupling is taken into account in this approximation. The Hamiltonian reads 



HmF — -ffload + -ffdot — ^MF 



\SM ) 



(36) 



where 



Jmf = J (1) + VJW2 + 3J2, 



(37) 



and 



tan <p = VsJ/ Jmf- 



(38) 



The last term in Hmf considers the restriction of Eq. (|l l|) . The expectation value of Hmf is minimized with respect 
to |S| 2 . The Kondo temperature can be estimated by 



T K = ttv\ J MF (H}|" 



(39) 
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using (S) determined by the self-consistent calculations (Appendix B). 

First let us consider the case in the absence of the Zeeman effect, E\m — -Es=i and Eoq = E$=i + A. The resonant 
level for the triplet state is threefold degenerate at Es=i = E$=i + A whereas the resonant level for the singlet state 
is at Eq = Eqq + A. These levels are separated by the energy A. The Lagrange multiplier A is determined to fulfill 
Eq. (|TT|) . Figure 4(a) shows the calculated results of Tk as a function of A. Both of Tk and A are in units of 
Do exp(— 1/j/Jmf)- We find that (i) Tk(A) reaches its maximum at A = 0, (ii) for A 3> Tk(0), Tk(A) obeys a power 
law 

T K (A) ■ A ta " 2 v — const., (40) 
and (iii) for A < 0, Tk decreases rapidly with increasing |A| and disappears at A = A c ~ — Tk(0); 

A c = -£>oexp(-l/W M F)(l + tan 2 ( / 3)(tanV) _sin2v - (41) 

These features are in agreement with the results of the scaling calculations. 

The behaviors of Tk(A) can be understood as follows. The inset of Fig. 4(a) schematically shows the Kondo 
resonant states. The resonance of the triplet state is denoted by solid lines whereas that of the singlet state is by 
dotted lines. (A) When A 3> Tk(0), the triplet resonance appears around /i whereas the singlet resonance is far above 
fi. (B) With a decrease in A, the two resonant states are more overlapped at fi, which raises Tk gradually. This 
results in a power law of Tk(A), Eq. (|40|). The largest overlap yields the maximum of Tk at A = 0. (C) When A < 0, 
the singlet and triplet resonances are located below and above \x, respectively, being sharper and farther from each 
other with increasing |A|. Finally the Kondo resonance disappears at A = A c . 

The conductance through the dot is given by 



G/(e 2 /h) = 



4F 1 F 1 / A 2 A 2 \ 4F 2 F 2 A 2 

^ L L L R / ^11 , ^10 \ , * L L 1 R ^00 



(ri + r^) 2 I (a - En) 2 + A 2 ! (s - E w f + A 2 ) (r 2 + Y\f (e - E 00 ) 2 + A 2 



(42) 



where T l a = TTv\V a _i\ 2 . The resonant widths are Ah/Aq = 2cos 2 tp/3, Aio/Aq = cos 2 tp/3, and Aoo/Ao = sin 2 ip with 
Ao = ttv\ Jmf(S) | 2 . The conductance G as a function of A is shown in Fig. 4(b), in a symmetric case of T L = T R 
(i = 1,2). G = 2e 2 /h for A > whereas G goes to zero suddenly for A < 0. Around A = 0, G is larger than the 
value in th£_.unitary limit, 2e 2 //i, which is attributable to nonuniversal contribution from the multichannel nature of 
our model.EJ 

In the presence of the Zeeman splitting, E\m — E$ = i — EzM, the resonant level of the triplet state is split into 
three. With increasing Ez, the Kondo effect is rapidly weaken except in the region of A ~ — E%. In Fig. 5(a), we 
show the Kondo temperature Tk in E%-A plane, in the case of ip = 0.157T. Figure 5(b) presents Tk as a function of 
A for several values of Ez- When Ez is large enough, the Kondo effect takes place only when the resonant state of 
|11) is overlapped with that of |00). Then Tk is the largest at A = — Ez and decreases with A being away from this 
value. At A = —Ez, Tk obeys a power law 

T K (A) • |A|l/(2+3tan 2 V ) = (43) 

which is indicated by a broken line in Fig. 5(b). This is qualitatively in agreement with the calculated results by the 
scaling method. 

Figure 6 indicates the conductance G in Ez-A plane, when ip — 0.157T and T l L = T l R (i = 1,2). G takes the value 
of 2e 2 /h around Ez — and A > 0, and also along the line of Ez = —A. (G > 2e 2 /h in the neighborhood of 
Ez = A = 0, as discussed above.) For sufficiently large Ez, our model is nearly equivalent to the anisotropic Kondo 
model with S = 1/2.E3'E3 Hence G = 2e 2 /h at A = -Ez and reduces to zero as A deviates from this value, in the 
same way as in Fig. 3(b) for the case of S = 1/2. 



V. CONCLUSIONS AND DISCUSSION 



The Kondo effect in quantum dots with an even number of electrons has been investigated theoretically. The Kondo 
temperature Tk has been calculated as a function of the energy difference A = Eqo — E$=i and the Zeeman splitting 
Ez, using the poor man's scaling method and mean field theory. The scaling calculations have indicated that the 
competition between the spin-triplet and singlet states significantly enhances the Kondo effect. When the Zeeman 
effect is irrelevant, Ez <S. Tk, Tk is maximal around A = and decreases with A obeying a power law. In a case of 
—A = Ez, the Kondo effect takes place from the degeneracy between two states, 1 00) and |11). Even in this case, the 
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contribution from the other states of higher energy, 1 10} and |1 — 1), plays an important role in the enhancement of 
Tk- As a result, Tk is maximal around Ez = and depends on Ez by a power law again. 

The mean field theory yields a clear cut view for the Kondo effect in quantum dots. Considering the spin couplings 
between the dot states and conduction electrons as a mean field, (fsM c k\)i we nri< ^ that tnc resonant states are 
created around the Fermi level fi. The resonant width is given by the Kondo temperature Tk. The unitary limit of 
the conductance, G ~ 2e 2 /h, can be easily understood in terms of the tunneling through these resonant states. In 
our model, the overlap between the resonant states of S = 1 and S = in the dot enhances the Kondo effect. The 
mean field calculations have led to a power law dependence of Tk on A and on Ez, in accordance with the scaling 
calculations. ■— . 

The mean field theory is not quantitatively accurate for the evaluation of Tk£9 (In the case of S = 1/2, the exact 
value of Tk is obtained accidentally.) In our model, the scaling calculations indicate that all the exchange couplings, 
jC 1 ), j( 2 ) ; and J, are renormalized altogether following Eq. ( |l6|) when |A| and Ez are much smaller than the energy 
scale D. In consequence two channels in the leads are coupled effectively for an increase in Tk- In the mean field 
calculations, the interchannel couplings are taken into account in Eq . (|37|) only partly. In fact, conduction electrons of 
channel 1 and 2 independently take part in the conductance, Eq. (|42|). By the perturbation calculations with respect 
to the exchange couplings, we find that mixing terms between the channels appear in the logarithmic corrections to 
the conductance.Eil We could improve the mean field calculations by adopting another form of the order parameter 
than Eq. (§|). 

Our calculated results explain the experimental findings by Sasaki et al.'XJ The Kondo effect is largely enhanced 
around A = when the Zeeman effect is irrelevant. The behavior of Tk in the presence of the Zeeman effect may 
be observed experimentally under higher magnetic fields. In experiments the value of the Zeeman splitting can be 
controlled by applying a magnetic field parallel to the quantum dot. More generally, it is possible to control several 
parameters in semiconductor quantum dots and to realize new situations which cannot be reached in traditional solid 
state context. The quantum dot systems, therefore, have the potential of tools to explore the Kondo physics further 
beyond the present theory. 
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APPENDIX A: MEAN FIELD CALCULATIONS FOR S = 1/2 

The original Hamiltonian for a quantum dot with one energy level reads 

h= J2 £ ^l M c aM + J2 Y,^ c lka d -+ R - c -)+ H ^ ( A1 ) 

a—L 7 R ka a—L,R kcr 

with 

H dot = 5^£b4<k + Ud^dl^. (A2) 

a 

For the state of one electron in the dot, the addition and extraction energies are given by E + = Eq + U — [i and E~ = 



fi — Eo, respectively. The parameters, Eq and U, in Eq. (A2) should be determined to fit these energies to experimental 



data. For conduction electrons in leads L, R, we perform a unitary transformation, c ka = (V£cL,ka + V R c R ,ka) /V , 
Cka = (-V R .iC L , krT + V L ,iC R ,kcr)/V, with V = ^J\V L \ 2 + \V R \ 2 , along the lines of Ref. §. We disregard the modes c ka 
which are uncoupled to the quantum dot. 

We consider the Coulomb blockade region for one electron, where both E + and E~ are much larger than the level 
broadening T = irvV 2 {y being the density of states in the leads^pand temperature. Integrating out the dot states 
with zero or two electrons by the Schrieffer- Wolff transformation,Lrtl we obtain the effective low-energy Hamiltonian 



H = ^e k c\ a c k(T + ^E a flf a + [s+4 u c feT + S-cl,^c kl + S z {c\,^c k] - c^'^ki) (A3) 

kk' 
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under a constraint of Eq. (p9|). In the second term we have included the Zeeman effect, Ef t ± — Eq ± Ez- The 
third term represents the exchange coupling between the dot spin and conduction electrons with J = V 2 / E c where 
1/E C = 1/E + + 1/E~ . By expressing the spin operator S as S+ — /j/j,, S- = fffi, S z = (ftfi - /j/j.)/2, one finds 
that Eq. (|Aj|) is identical to Eq. ((28^. 

The mean field Hamiltonian, Eq. (|3l|), includes "energy levels" for pseudo-fermions, E a = E a + X, which are coupled 
to the leads with "tunneling amplitude," V — — \/2J(S). The Green function for the pseudo-fermions is 



G a (e) 



e - E a + iA ' 



(A4) 



where A = 7w|U| 2 . This represents the resonant tunneling with the resonant width A. 
The expectation value of the Hamiltonian, Eq. (|3l|) , is written as 



E 



MF 



E 



A E„ 



tan 



E„ 



A £2 + ^2 

— In s 

2tt L> 2 



■KVJ 



(A5) 



where Z?o is the bandwidth in the leads.0 We set fi = in this appendix. The constraint of Eq. (29) is equivalent to 
the condition 



OEmf 



- Vtan" 1 -^--1 = 0. 

IT J? 



dX 7T ^ E a 

This yields Eq + A = 0. The minimization of Emf with respect to A (or |(S)| 2 ) determines A 

E 2 +A 2 



OEmf 
dA 



^E^ 



n 2 



1 



Q. 



For Ez =0, we find 



A = D exp[-l/2i/J] = A . 



This is equal to the Kondo temperature, Tk- For ^ 0, Eq. (|A7| ) yields 

A 2 +E% = Aq. 

Using the T-matrix, T, the conductance through the dot, G, is given by 

,2 



G=j(2ttv) 2 J2 \{R,k'a\f\L,ka)f 



e%g Tr ,2 \Vl\ 2 \Vr\ 2 



4r L r fl . 



E 



cr 

A 2 



h (T L + T R ) 2 ^ ( e - ^ ff )2 + a 2 



(A6) 



(A7) 



(A8) 



(A9) 



(A10) 



where T a = nv\V a \ 2 . This yields Eq. (||) in the text. On the second line in Eq. jAlOj), |i/> fe(T ) = cj. CT |0) = (Vl|Z, ka) 
Vr\R, ka))/V, and the T-matrix is evaluated in terms of the Green function, Eq. (A4), |U| 2 G<j(£ = £&). 



APPENDIX B: MEAN FIELD CALCULATIONS IN THE PRESENT MODEL 



For the spin states of the coupling between the spin triplet S = 1 in the dot and a conduction electron, we introduce 
spinors of £ = 1/2 and 3/2. Using the Clebsch-Gordan coefficients, they are given by 
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^ 2 = E 



n 



(0 

3/2 



E 



(^cg + AVg)/^ \ 

(^/ 1 t _ 1 cg+/ 1 t oC g)/^y 

' Jii c fc| _ > 

(-AVS + ^AVg)/^ 
(/t_ lC « - V5/i,c«)/V5 



V 



(Bl) 



(B2) 



The exchange couplings between the triplet state and conduction electrons, Eq. (U3), can be rewritten as 



_9(=)('0to(0 i rjWtoW 
Z "l/2"l/2 T "3/2 3/2 



(B3) 



In the same way we define the spinors of S = 1/2 to represent the spin couplings between the singlet state 5 = and 
a conduction electron 



f t m 

■>00 c fc|._ 

_ f \ (0 



(B4) 



where i = 2 and 1 for i — 1 and 2, respectively. The conversion between the triplet and singlet states, Eq. ([13]), is 
rewritten as 



(B5) 



In H s_1 + ^j 5 '- 1 ^-*0 j a mode of the largest coupling with S = 1/2 is given by 

3 = cos ip(l^ 2 + sin ^vjK 1 ) 



(B6) 



for > j( 2 ), which is Eq. (|35| ) in the text. The corresponding eigenvalue is given by Eq. (|37| ) and </? is determined 
as in Eq. (|38|). 

The mean field Hamiltonian, Eq. (^6|), represents the resonant tunneling through the energy levels for the pseudo- 
fermions, Esm — Esm + A. The expectation value of Eq. (^6|), £mf> is evaluated in the same way as in Appendix A. 
dE MF /d\ = yields 



, -i An . , _i Aio _! A o 
tan — h tan — h tan — — = 7r, 



(B7) 



En Eiq Eqq 

where the resonant widths are An/Ao = 2cos 2 < / 9/3, Aio/Aq = cos 2 </?/3, and Aoo/Ao = sin 2 99 with Ao = 
~kv\ Jmf(3) | 2 . We set /i = here. Minimizing Emf with respect to Ao, we obtain 



I cos 2 <p In g " + A " + - cos 2 p In " '" 
3 £>o 3 



p2 , A 2 p2 1 A 2 

[i ' i° + sin Vn 00+ 00 



D 2 



n 2 



Equations ( |B7| ) and ( |Bq ) determine A and A (or |(S)| 2 ). 
The conductance through the dot is given by 



(B8) 



G= -{2-Kvf 
h 



J2 \{R,k'a',j\f\L,ka,i 

r 



E 



1 g 



= °-(2^) 2 . 



L |T<I„/,«\|2 



(B9) 



where Tl, = ixv \V a .tf and = (V L ,i\L,ka,i) + V R ,i\R,ka,i))/Vi. The T-matrix can be evaluated, using the Green 

function for the pseudo-fermions, Gsm(s) = [s— Esm + iAsm] -1 ) a s in Appendix A. This yields Eq. (E3) in the text. 
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Figure captions 

Fig. 1: (a) The energy diagram for the spin states \SM) considered in our model. A = E o — Es=i and Ez is 
the Zeeman splitting, (b) Spin flip processes between the spin states. The exchange couplings jW involving the 
spin-triplet state only are accompanied by the scattering of conduction electrons of channel i. Those involving both 
the spin-triplet and singlet states, J, are accompanied by the interchannel scattering of conduction electrons. 

Fig. 2: The scaling calculations of the Kondo temperature Tk as a function of A, (a) when the Zeeman splitting 
is irrelevant, Ez <C Ik, and (b) in a case of Ez = —A. D is the bandwidth in the leads. In both the figures, a, 
6/ir = 0.25; b, 0.15; and c, 0.10 where tan 2 9 = jM/jW. 

Fig. 3: The mean field calculations for the Kondo effect in a quantum dot with S = 1/2. (a) The Kondo temperature 
Tk and (b) conductance through the dot, G, as functions of the Zeeman splitting Ez- 7k and Ez are in units of 
D exp(— l/2i/J) and G is in units of (2e 2 /h) ■ 4:TlTbJ(Tl + T^) 2 . Inset in (a): The Kondo resonant states created 
around the Fermi level fi in the leads, (A) in the absence and (B) presence of the Zeeman splitting. The resonant 
width is given by Tk- 

Fig. 4: The mean field calculations for the Kondo effect in the present model. The Zeeman splitting is disregarded 
(Ez <C Tk). (a) The Kondo temperature Tk and (b) conductance through the dot, G, as functions of A = Eqq —Es=i- 
Tk and A are in units of D exp(— 1/v Jmf)- G, in units of 2e 2 /h, is evaluated in a symmetric case of T l L = T l R (i = 1, 2). 
tanip = V3J/Jmf where a, ip/ir = 0.25; b, 0.15; and c, 0.10. Note that ip/n < 1/6 in this approximation (case a is 
only for reference). Inset in (a): The Kondo resonant states for 5 = 1 (solid line) and for S = (dotted line) when 
(A) A > T K (0), (B) A - T K (0), and (C) A < 0. 

Fig. 5: The mean field calculations for the Kondo temperature Tk in the present model, (p/n = 0.15 where 
tan (p — \/3J/Jmf- All of Ez, A, and Tk arc in units of D cxp(— 1/^Jmf)- (a) Tk is plotted in Ez-A. plane, by 
contour lines drawn every 0.25. The brighter color indicates the larger values of Tk- (b) Tk as a function of A when 
E z is fixed at a, 0; b, l; c, 2; d, 5; and e, 10. The broken line indicates T K in the case of — A = E Z - 

Fig. 6: The mean field calculations for the conductance G in the present model. G is plotted in Ez-A plane, by 
contour lines drawn every 0.2 x (2e 2 /h). The brighter color indicates the larger values of G. Ez and A are in units 
of D exp(— 1/vJmf)- vl™ = 0-15 where tany = V3J/Jmf, and T l L = T Z R (i = 1,2). 



13 



(a) energy 



(b) 



■oo 



: S=1 



A 



|00> 

I1-1> 

|10> 
|11> 




Figure 1 (M. Eto and Yu. V. Nazarov) 
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